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1 . Introduction 

With  the  ability  to  solve  general  nonlinear  optimization  (programming) 
problems  has  come  the  development  of  more  sophisticated  and  realistic  models. 

One  example  of  this  in  the  logistics  area  is  the  presentation  in  [1]  of  an 
optimal  inventory  model.  This  model  was  solved  by  the  Sequential  Unconstrained 
Minimization  Technique  (SUMT)  using  as  the  interface  between  the  algebraic 
statement  of  the  problem  and  the  canned  algorithm,  four  FORTRAN  subroutines. 

This  traditional  mode  of  input  to  the  nonlinear  programming  algorithm  was 
required  since  function  value  and  derivative  information  is  necessary  to 
solve  the  problem  efficiently.  A first  attempt  at  providing  a simpler  model/ 
algorithm  interface  was  done  by  Pugh,  e.t  al.,  [2]  and  was  a very  sophisticated 
computer  code  allowing  for  general  forms  of  symbolic  input.  Unfortunately, 
this  code  was  slow.  A more  sparse  approach  to  problem  representation  was 
developed  in  [3].  Internally,  the  code  was  fast,  but  the  user  was  restricted 
to  numerical  input.  This  code  Is  in  the  process  of  undergoing  revision  [A] 
and  in  this  memorandum  the  inventory  model  given  in  [1]  is  presented  in  the 
new  format  which  allows  for  mnemonic  symbolic  input  for  variable,  constant, 
and  function  names.  This  will  make  the  existing  nonlinear  programming  codes 
more  accessible  to  practitioners.  This  memorandum  is  an  adjunct  to  the  users 
manual  [A]  and  should  be  used  in  conjunction  with  it  and  the  original  paper  [1]. 
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2 . The  Inventory  Model 

Below  is  given  a brief  development  of  the  optimization  problem. 

For  more  details  the  reader  is  referred  to  (IJ  and  [5]. 

Fixed  Inputs  and  Constants 

N = number  of  items  in  inventory  (n  = 2N  is  the  number 

of  variables  or  unknowns  in  the  nonlinear  programming 
problem) 

c.  -■  item  unit  cost  (in  dollars) 

1 

V.  = item  unit  volume 
1 

= mean  demand  per  unit  time 

(p.,0.)  = mean  and  standard  deviation  of  lead  time  demand 
distribution  for  1th  item. 


Let 

the  normal  density  function,  and 

the  right  hand  cumulative  normal 
distrlb\ition  function. 

= investment  limit  in  dollars 

K2  = reorder  work  load  limit 

= volume  limit 


1(1  (x)  = ( 1 / e 
<i(z)  = <i(x)dx 


Derived  Quantities 

= \lol  + (z-p.)^] 


0 . 
1 


2 


(z-p . ) 
1 


EOH.  = r.  + 


Qj/2  - 


h.  + B. 


(Qi.ri) 


2 


i- J.-)/ 


The  optimization  problem  is 


,N 


subject  to 


unimize  ^ , B,(Q.,r.)  (expected  time-weighted  shortages) 

(Qi-ri)  i 1 1 


|r.  + ■^  ” *^i  — *^1  (Investment  constraint) 

11  ^i 

^ ^ (reorder  workload  constraint) 

|r^  ^ ^ ~ ^'i  ^ *^i^*^i’*^i^|  — *^3  ’ (capacity  constraint) 


> 0 , i = 1 N 


3 . The  Model  in  Factorable  Form 

The  factorable  programming  language  which  provides  the  interface 
between  the  algebraic  statement  of  the  optimization  problem  and  canned 
optimization  algorithms  assumes  that  the  problem  can  be  written  in  the 
following  canonical  factorable  form: 

min  X (x) 

X C e" 

Subject  to:  ^ X^(x)  < , for  i=l,...,N-l 

(possibly  = -»>  and/or  = +^)  where  X^(x)  i x^  (for  1 = l,...,n) 

and  the  remainder  are  defined  recursively  as  follows: 

given  X*’(x),  for  p = l,...,i-l,  then  for  i = n+l,...,N  , 

. i-1  i-1  p . 

X (x)  = I TMxP(x)]  + I r.  [xP(x)]  • (x‘’(x)]  , 

p=l  P p=l  q=l  P’'^ 

where  the  T's,  U's,  and  Vs  are  functions  of  a single  variable. 

- 3 - 
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In  writing  down  his  particular  problem  in  canonical  factorable 

form  the  user  should  obey  two  fundamental  rules.  Tlie  quantities  X^(x) 
are  called  concomitant  variable  functions  (cvfs)  and  on  the  right-hand 

side  of  the  defining  equation  for  any  there  cannot  appear  any  X^ 

where  p ^ i . Furthermore,  no  more  than  two  factors  can  appear  in  any 
product  term. 

The  new  format  allowed  for  the  factorable  programming  language 
allows  for  symbolic  variable  names  Instead  of  just  numbers  for  the  cvfs. 
The  user  may  use  any  mnemonic  devices  to  help  in  setting  up  the  factorable 
canonical  form.  For  the  problem  example  above  this  could  take  the  follow- 
ing form: 


z . = 
1 

(ri-Pi) 

, i = 1,  . 

..,N 

(r.-p.) 

II 

Z.  + Q. 
1 1 

, i = 1,. 

..,N 

(ri^l-U^) 

U . = 

1 

Zf/Oi 

, i = 1, . 

..  .,N 

(ri-pJ/Oi 

''i  = 

Vl/Oi 

, i = 1,. 

..,N 

(r.4<j.-p.)/a. 

Wi  = 

2 ^ 2 
0 . + z . 

1 1 

, i = 1 , . 

..,N 

(cJ+(ri-Uj)^ 

2 ^ 2 
Of  +y, 

. 1=1,- 

..,N 

(o^+(r.-H}.-u.)2) 
1 111 

, u.  2,» 
p.  = — / ^ e"*^  '^^dt 

" /?■  ™ 


1 r^i  -t  /2,^ 

q = J „ e dt 


Di(ri,Qi>  = 1/2  w^H-p^l  - z.<(>(u^) 

1 '’i 

- j t^[l-q^]  + ~ y^4»(vj),  i = 1, 


iP(r^)  - 
3j(r^+Q^) 
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(Q^.r.)  = 

D.(r 

,g.)/Q.  , 

i - 1 N 

1 1 

1 

L ^ 1 

EOH^  = 

r.  + 
1 

Qi/2  - + R^(Q.,rJ, 

i = 1,...,N  (EOH^) 

EXPINV  = 

Z 

^i=l 

c.  EOH. 
1 1 

(EXPECTED  INVESTMENT) 

EXPVOL  = 

z '' 
1=1 

V.  EOH. 
1 1 

(EXPECTED  VOLUME) 

REORD  = 

z "" 
^i=l 

X./Q. 
1 1 

(WORK  LOAD) 

SHORT  = 

B.(Q.,r.) 

(EXPECTED  TIME- 

1=1 

1 1 1 

WEIGHTED  SHORTAGES) 

Under  the  new  system  the  user  first  declares  the  names  of  his 
variah’  and  then  the  names  of  the  data  constants  he  uses.  He  also  has 
the  lity  of  placing  bounds  on  ranges  of  variables  using  their  symbolic 

n 

In  Appendix  A is  given  the  listing  of  the  card  input  for  the  inventory 
model  developed  in  [1].  A slightly  more  elaborate  model  was  presented  in  [5] 
based  on  material  in  [1].  The  card  listing  for  this  model  is  given  in 
Appendix  B. 

A . Conclusion 

A brief  description  of  the  manner  in  which  a particular  logistics 
optimal  Inventory  model  can  he  described  in  the  new  factorable  programming 
interface  has  been  given.  Hopefully  the  new  flexible  format  which  allows 
for  symbolic  names  of  variables,  constants,  and  functions  will  make  the  use 
of  optimization  codes  more  widespread. 
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APPENDIX  A 


SYMBOLIC  FACTORABLE  PROGRAMMING  INPUT  FOR 
SIMPLE  SCHRADY-CHOE  MODEL  DESCRIBED  IN  [I 
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APPENDIX  B 


SYMBOLIC  FACTORABLE  PROGRAMMING  INPUT  FOR  MODIFIED 
SCHRADY-CHOE  MODEL  DESCRIBED  IN  [5] 
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